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The possibility that the apparent room temperature ferromagnetism, often measured in Co-doped 
ZnO, is due to uncompensated spins at the surface of wurtzite CoO nanoclusters is investigated by 
means of a combination of density functional theory and Monte Carlo simulations. We find that 
the critical temperature extracted from the specific heat systematically drops as the cluster size is 
reduced, regardless of the particular cluster shape. Furthermore the presence of defects, in the form 
of missing magnetic sites, further reduces Tc- This suggests that even a spinodal decomposed phase 
is unlikely to sustain room temperature ferromagnetism in ZnOiCo. 



In recent years the search for ferromagnetism in 
insulating oxides doped with small quantities of 
transition metals has become a topic generating much 
debate in the literature. Taking ZnO: Co as the 
proto- typical example for this class of materials, many 
experimental groups have reported room temperature 
ferromagnetism^ ^ , whereas several other have failed to 
find any such evidencd^. Notably, growth conditions, 
sample morphology and spatial Co distribution play a 
crucial role in determining the magnetic properties. In 
particular, there is now an emerging view that samples 
with high structural quality and uniform Co distribution 
do not result in long-range ferromagnetism at high 
temperature, and that the magnetism may be related to 
structural^ or point defects^. 

Given the unsettled experimental landscape it should 
not surprise that a number of interesting and competing 
theoretical models have been proposed. In general 
explanations involving standard mechanisms for the 
magnetic interaction appear problematic. Schemes 
leading to short range magnetic coupling such as super- 
exchange need a Co concentration, [Co], exceeding the 
percolation threshold. This is around 20 % for an 
interaction extending to the nearest neighbour sites of 
an fee lattice, and it is about a factor of 5 larger than the 
typical experimental concentrations. Similarly carrier 
mediated mechanisms are not sustained by experimental 
evidence, which shows both paramagnetism in presence 
of abundant free carrierS^^and ferromagnetism deep in the 
insulating region of the phase diagrarrP. More generally 
carrier concentration and mobility have an apparent little 
correlation to the magnetic propertied. Furthermore 
carrier mediated mechanisms are difficult to validate 
on a solid theoretical ground by using first principles 
calculations, since the empty Co d levels are usually 
erroneously predicted too shallow at the edge of the ZnO 
conduction band^. 

Thus one has to look for more complex mechanisms 
for the magnetic interaction. Amongst these, the donor 
impurity band exchange model (DIBE)^*^ has enjoyed 
considerable popularity in the experimental community. 
According to the DIBE the magnetic interaction among 
Co^+ ions is mediated by donors, whose charge density 
is localized over large hydrogenic orbitals, so that 
the relevant percolation threshold becomes that of the 



donors and not that of the Co. Unfortunately the 
model still fails at the quantitative ground, since room 
temperature ferromagnetism needs prohibitively large 
exchange coupling between the donors and the Co^" 
A sec ond, recently proposed scheme, is the two-species 
modePESI in which Co-oxygen vacancy pairs (CoV) 
act as a second magnetic center in addition to Co^+. 
Interestingly CoV can interact magnetically strongly up 
to the third nearest neighbours, where the percolation 
threshold drops to 7 %. Although this is still too high 
to justify a long-range ferromagnetic order, it suggests 
that room temperature ferromagnetism can be achieved 
in samples where Co ions are non-uniformly distributed. 

The extreme limit of Co segregation is represented 
by the formation of some secondary Co-based phase. 
Usually metallic Co is excluded by X-ray dataP^. 
Unfortunately all other compounds in the Zn-Co-0 
phase diagram (CoO, C02O3, C03O4, ZnCo204) are 
either non-magnetic or antiferromagnetic with low Neel 
temperatures, so that there is no obvious candidate 
material to form ferromagnetic clusters in ZnO. Still, 
recently Dietl et al. suggested that the apparent 
room temperature ferromagnetism in ZnO:Co may be 
due to uncompensated spins at the surface of CoO 
nanoclusters maintaing the ZnO wurtzite (WZ) lattice 
structureP. Indeed both W Z and z incblende (ZB) CoO 
were synthesized in the bullP^E^E^, but little magnetic 
characterization exists. Unfortunately calculations based 
on density functional theory (DFT) and Monte Carlo 
(MC) methods are discouragin^^^. In fact all the possible 
CoO polymorphs, including WZ and ZB, do not display 
a ferromagnetic order but instead they are characterized 
by a general spin frustration. As a result, the critical 
temperatures, Tc, extracted from the specific heat are 
well below room temperature. Notably such a frustration 
was recently confirmed experimentally^^ for the WZ 
phase. 

However, even if one can exclude magnetism in the 
bulk, the question of whether or not uncompensated 
spins can order at the surface of a nanocluster remains 
open. In particular this is an intriguing question since 
the dominant exchange parameter extracted from DFT 
for wurtzite CoO, the one giving frustration, is large. 
One can then speculate that reducing frustration (as 
it happens on a surface) can help in enhancing Tc- 
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Such an open question is investigated here, where we 
use MC simulations for extracting the thermodynamical 
properties of WZ CoO nanoparticles of different shapes 
and dimensions. 



I. FOUNDATION OF THE MODEL AND 
COMPUTATIONAL DETAILS 

We perform MC calculations for a classical Heisenberg 
Hamiltonian of the form 

F=-i^J,,S^.5;+^i?(s..n)' , (1) 

i,j i 

where Si is a classical spin located at site i {\Si\ = 3/2 
for Co^"*") and Jij is the exchange parameter between 
spins at sites i and j. The second term in ([ij describes 
the hard-axis easy-plane crystal anisotropy, n is a unit 
vector along the WZ c-axis and the zero-field splitting 
for Co^+ taken from electron paramagnetic resonance'^- 
is D = 2.76 cm-i. 

The values for the various exchange parameters, J, for 
bulk WZ CoO have been calculated previously' from 
DFT using the LDA-I-J7 extension of the local density 
approximation (LDA) . The procedure used was to fit the 
total energies of a number of reference DFT calculations 
for different magnetic supercells to the Hamiltonian of 
equation ([T]). For WZ CoO four different J's are enough 
to describe the dominant magnetic interaction, with an 
error over the reference DFT calculations smaller than 
1 meV/Co. We have then found: Ji = 6.1 meV, J2 — 
—36.7 meV, J3 — —0.2 meV, J4 — —5.2 meV, where 
the index n in J„ refers to the neighbour degree (i.e. 
"1" is for first near neighbours Co^+). The dominant 
interaction, J2, is thus for first nearest neighbours in the 
{001} plane (second nearest neighbours overall) and it is 
antiferromagnetic. This is the origin of the frustration 
and of Tc's below room temperature for the bulk. The 
same parameters are used for the calculations of finite 
particles presented here. 

The ground state of a given particle is found by the 
simulated annealing method, since frustration introduces 
many low energy configurations differing considerably 
from the ground state but with minor energy differences 
from it. Such an energy landscape clearly makes the 
conjugate gradient and the steepest descent schemes 
ineffective. The temperature dependence is then 
investigated with Monte Carlo simulations. We use 
the standard Metropolis algorithm where the acceptance 
probability of a new state is 100 %, if the new state 
has an energy lower than that of the old one, and it 
is given by the Boltzmann factor otherwise. For each 
different magnetic cluster we equilibrate the system at 
a given temperature and then extract thermodynamical 
quantities by sampling over several millions MC steps. 
In particular we extract Tq from the peak in the specific 
heat, C. This becomes necessary since an obvious order 



parameter is difficult to find for these highly frustrated 
clusters. 

Such a computational scheme is not free of uncertainty. 
Firstly, the reference DFT calculations are dependent on 
the specific choice of exchange and correlation function 
used. LDA+U is certainly suitable for CoO and our 
parameterization of the U — J parameter is based on 
total energy considerations, i.e. on fitting the structural 
properties and not the band-structur#3. Secondly the 
Heisenberg model used contains only J's extending over a 
limited range and includes only pairwise interaction. The 
first approximation appears acceptable given the good 
quality of the fit to DFT, while the second one is more 
difficult to assess. Nevertheless we have used rocksalt 
CoO to estimate the error and find that typically our 
Tc's for the bulk are underestimated by about 30 /^f^. 
Considering that the Co valence in rocksalt and WZ CoO 
is the same, we speculate that the same error found for 
the rocksalt phase can be transfered to WZ. 

Still the uncertainty of using bulk parameters for finite 
clusters simulations remains. However, it is important 
to bare in mind that here we do not consider free- 
standing CoO nanoparticles, but instead CoO clusters 
embedded into a ZnO matrix. This means that the local 
chemical coordination of each Co atom (i.e. the fourfold 
coordination to O) in the cluster is identical to that of 
bulk CoO. Thus the finite size affects only the magnetic 
coordination and one can safely use bulk parameters for 
clusters. This, of course, would be inadequate in the 
case of free surfaces or grain boundaries for which a new 
parameterization is needed. 



II. RESULTS 

If present, it is likely that CoO nanoclusters form in 
ZnO either during the growth process or during post- 
growth treatment. DFT calculations in fact confirm 
that there is an energy gain when moving Co ions to 
nearest neighbouring positionJ^ so that clustering is 
highly probable. This can possibly be tuned by growth 
parameters and eventually the presence of additional 
dopants^. Still, even assuming that ZnO:Co during the 
growth has enough kinetic energy to form large CoO 
clusters, it is hard to predict whether such clusters should 
have a preferential shape. Indeed simulations of spinodal 
decompositiorP2l seem to suggest the formation of 
highly Co-rich regions without any particular geometrical 
structure. This is then a different situation from that of 
free-standing nanoparticles, where the presence of free 
surfaces drives the particle geometry. For this reason 
we have looked at three different particle shapes, namely 
spherical, cylindrical and perforated spherical. The last 
ones are spherical particles where a fraction of the Co 
sites is randomly removed. 
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1. Spherical Nanoparticles 

We consider spherical particles first. These are 
constructed by simply removing all the atoms beyond 
a sphere of a give radius, R. In doing so the position 
of the center of the sphere influences the atomic details 
of the external surface so that different particles with 
the same radius can be made. We flnd that all the 
particle properties are relatively sensitive to the position 
of the particle center (chosen within the WZ unit cell) 
for small particles, but they become progressively center- 
independent as the radius gets larger. This is expected 
since larger surfaces allow all the thermodynamical 
quantities to self-average. For this reason we perform 
averages over the position of the center only for the 
smaller particles investigated. 

i"/^-,^ \"^^s- 

^-'/W \-/\^/ 1-/ W 
~ /l"V^-,*'\--/V/^-/V 
FIG. 1: [Color on line] Ground state spin configuration of 
spherical particles of different size: 12 A (top) and 18 A 
(bottom). Note the large degree of compensation in the inner 
core of the particles and the presence of non-compensated 
spins at their surfaces. 

We start our analysis by presenting the ground 
state spin configuration as calculated with simulated 
annealing. This is shown in figure [T] respectively for a 
small (R — 12 A) and a large (i? = 18 A) particle. 
In the figure the particles are oriented with the WZ c- 
axis pointing perpendicularly to the page, so that the 
a-b planes are visualized. The figure displays a clear 
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TABLE I: Table listing the uncompensated magnetic moment 
per Co in the ground state, /i, and the critical temperature, 
Tc, extracted from the specific heat, for spherical particles 
of different radius. We also report the number of Co atoms 
contained in the particle, A^'co, and both particle volume, V , 
and area of the surface, S. 



spin compensation of the particle inner core. As for the 
case of bulk WZ CoO, the strong first nearest neighbour 
antiferromagnetic constant in the a-b plane, J2, drives the 
frustration. As a result the spins in the a-b plane align at 
120° with respect to each other and the net moment per 
plane vanishes. At the surface the frustration is lifted by 
translation symmetry breaking and a net uncompensated 
moment, fi, emerges. Its direction and intensity depends 
on the particle size, and for small particles on the details 
of the surface geometry. In any case the moment is always 
rather small and it is then difficult to visualize by simply 
looking at figure [T] 

In table |T] we list such a ground state uncompensated 
magnetic moment (per Co atom) for spherical particles 
of different sizes. In the table for each radius we report 
the upper value obtained over a number of different 
geometrical realizations of the particle. In general 
the uncompensated magnetic moments are rather small, 
so that the cluster model itself is at best capable of 
explaining only weak magnetism in ZnO:Co. The table 
also allows us to extract an upper limit for the magnetic 
moment at room temperature due to WZ CoO particles. 
In fact, by assuming that magnetization reversal is driven 
by coherent rotation we estimate that particles with 
radii larger than 15 A (containing about 800 Co atoms) 
are large enough to be superparamagnetically blocked 
at room temperature. Their ground state (at T = 0) 
magnetic moment is calculated between 0.1 (Ib/Co and 
0.05 /xb/Co (Table |l]). Thus, if magnetism originates 
from CoO clusters with uncompensated spins, then a 
ferromagnetic signal at room temperature cannot be 
associated to magnetic moments in excess of 0.1 /ib/Co. 
In fact particles with larger n have a radius smaller than 
15 A, they are not be superparamagnetically blocked at 
room temperature, and hence they cannot contribute 
to the ferromagnetic signal. However, even for the 
possibility of weak ferromagnetism to be sustained one 
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has to demonstrate that the small magnetic moments 
survive at room temperature. 

Unfortunately for these frustrated system the magne- 
tization is not a good order parameter and no sublattices 
magnetizations can be identified. Moreover both the 
uncompensated moment and the particle susceptibility 
turn out to be rather noisy quantities, so that no Tq 
can be extracted from them. We then calculate Tq from 
the analysis of the specific heat, C, as a function of 
temperature. This is shown in figure [2] for a number 
of nanoparticles of different radii. The upper curve 
corresponds to i? = 6 A and the lower to i? = 
30 A with the curves in between corresponding to radii 
incrementing by 1 A. From the curves we can clearly 
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FIG. 2: The specific heat, C, as a function of temperature for 
finite clusters of wurtzite CoO. The curves are for particles of 
difi'erent radii, ranging from 6 A (top) to 30 A (bottom), in 
increment of 1 A 

identify a peak that becomes more diffuse as the particle 
size grows. Nevertheless we can still assign to the peak 
position the Tq of each particle. These are also reported 
in Tab. |l] Interestingly the Tq calculated for bulk WZ 
CoO with the same parameterization used here is only 
160 K^^, i.e. it corresponds to a particle with a radius of 
16 A. 

Finite size scaling theo ry for the Heisenberg model and 
ferromagnetic interactior r^ I predicts a variation of the 
critical temperature with particle radius of the form 

T^-Tc{R) RV^' ^ (2) 
Too \RoJ 

where _Ro and Too are respectively the correlation radius 
and the critical temperature for the bulk. This cannot 
be applied directly to our results since Too is smaller 
than Tc{R) for large clusters, i.e. the left-hand side of 
equation (|2| becomes negative for some R. We have 
then fitted the calculated Tc{R) to the equation (|2|, 
by assuming a different T^o- The results are reported 
in figure |3] In the curve we show two fits obtained by 
taking Too = 200 K and performing the fit respectively 
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FIG. 3: [Color on line] Critical temperature as a function of 
the particle radius R (* symbols). The lines are fits to finite 
size scaling theory [Eq. ^]. The solid red line corresponds 
to the fit obtained for i? > 6 A , while the dashed blue curve 
is for i? > 10 A. 



over the range i? > 6 A (solid red curve) and i? > 10 A 
(dashed blue curve). As expected the fit improves when 
particles with small radii are excluded, since the scaling 
low is strictly valid in the vicinity of Too • In any case we 
clearly observe that a rather good fit can be obtained for 
critical exponents in the range of that expected for the 
Heisenberg model (/3 ~ 1.42), when Too is about 200 K. 
For these values we find a correlation radius of the order 
of 6 A. 

We are at this time uncertain of why Too calculated 
from MC for an infinite system (in practice it is 
calculated for a finite system with periodic boundary 
conditions and sufficiently large cells) differs from that 
extrapolated from finite size scaling. In general C(T) 
for finite particles is much more diffuse than that 
calculated by using periodic boundary conditions, so 
that the assignment of Tq may be affected by some 
errors. However, such an uncertainty is smaller than 
the difference between the two temperatures (~40 K), 
so that an alternative explanation is needed. We believe 
that the surface and the core of the particles contribute 
in a substantial different way to the Tq so that scaling 
theory cannot be applied directly. In fact, depending on 
the magnetic coupling at the surface compared to that 
of the bulk, it was already demonstrated that Tq can 
eit her incr ease or decrease as a function of the particle 
gj^^HHlllll T]-^jg nieans that depending on the details of 
the magnetic interaction T^ as a function of the particle 
size can approach Too asymptotically either from below or 
from above. In the present case of frustrated interaction 
it appears that the bulk value is approached in a non- 
monotonic fashion so that there exits spherical particles 
with a Tq larger than that of the bulk. This aspect will 
be further investigated in the next section. 

In any case we find that spherical particles, despite 
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FIG. 4: [Color on line]. Specific heat as a function 
of temperature for a representative sample of cylindrical 
nanoparticles with different radii and lengths. 



for large radii present a Tq larger than that of the 
bulk, still fall quite short from being magnetic at room 
temperature. Furthermore the particles with the largest 
magnetic moment are those showing the lowest Tc, 
which makes us concluding that spherical wurtzite Co 
particles cannot be at the origin of the room temperature 
magnetism in ZnO:Co. 



length, he- For instance he is 2 A, 5 A, 10 A and 20 A 
respectively for R being 2 A, 5 A, 10 A and 15 A. Thus it 
is tempting to propose the empirical relation ~ R- As 
a consequence for particles with R/h < 1 there is only one 
relevant dimension, R. In fact for long cylinders {R <C h) 
the surface to volume ratio remains constant at 2/R, so 
that the magnetic energy density does not change with 
the cylinder length. This seems to be at the origin of the 
saturation of Tc with h. 

Interestingly for R = h the area of the surface of the 
cylinder and that of the sphere become identical, while 
the volumes follow the relation T4phcrc — 4/3V^yiindcr- 
In general we expect that the Tc of a nanoparticle has 
both a surface and a volume contribution. In the case 
R — h the surface contribution is identical for spheres and 
cylinders, while the volumetric contribution is expected 
to be larger for the cylinders, since these have a larger 
volume. Thus we expect that the Tc of a sphere of radius 
R is lower than that of a cylinder of radius R and h — R, 
as indeed confirmed by comparing the tables I and II. 

Finally we take a look at the magnetization, finding 
that this is always small, typically < lO^'^/ie/Co, and 
varies little with either changes in the particle dimensions 
or with the temperature. In summary, from our 
analysis it does not appear that cylindrical nanoparticles, 
similarly the spherical ones, are supportive of room 
temperature ferromagnetism. 



2. Cylindrical Nanoparticles 

In order to further investigate the relative importance 
of the surface and bulk contributions to the magnetism 
we consider cylindrical nanoclusters, constructed with 
the cylinder axis oriented along the WZ c-axis. By 
varying both the length and radius of such a particle 
the surface area may be considerably altered whilst the 
total volume and hence the total number of atoms in 
the cluster remains constant. In figure |4] we present 
the specific heat for a number of cylinders of different 
dimensions, where we can still clearly observe the 
presence of a peak. As in the case of spherical particles we 
associated the critical temperature to the peak position. 

These are presented next in table II, where we list also 
the particle volume and the area of the surface. The 
most relevant feature emerging from the table is the 
rather sharp dependance of Tc on the particle radius 
and its insensitivity over its length. This results in 
the interesting finding that particles presenting the same 
volume but different aspect ratios can display rather 
different Tc's. For example the volume of a particle with 
R = 5 A and /i = 30 A is about 30% larger than that of 
a particle with i? = 10 A and ft, = 5 A. Nonetheless its 
Tc is a factor 3 smaller (44 K against 122 K). 

In addition we find that the critical temperature of 
cylindrical nanoparticles becomes almost independent 
from the length of the cylinder beyond a certain critical 



3. Perforated Spherical Nanoparticles 

Since empty Co sites in a CoO nanoparticle may serve 
to lift the spin frustration by eliminating a fraction of 
the neighbours of each Co atom, it is worth analyzing 
the dependence of Tc over the Co concentration. This 
essentially corresponds to studying a random ZnO:Co 
alloy with large Co doping. In general we expect 
that, as long as the Co concentration exceeds the 
percolation threshold for the wurtzite lattice (19 % for 
nearest neighbour interaction), a magnetic order will 
be found. The question remaining is whether or not 
magnetization and Tc will increase with respect to their 
values for stoichiometric CoO when the Co concentration 
is reduced. 

The nanoparticles investigated here are constructed by 
taking one of the previously made spherical nanocrystals 
and then removing a chosen fraction of magnetic ions 
at randomly chosen sites. We consider particles with 
a radius of 19 A, which are large enough to show a 
substantial Tc, but not enough for the thermodynamical 
properties to saturate at their bulk values, i.e. the surface 
of the particles still makes a non-negligible contribution. 
Furthermore we find that particles with smaller radii 
are too sensitive to Co vacancies, to the extent that in 
general at Co concentrations below 90% a Tc is not easily 
located. Note that, although we name Co vacancy a site 
where the Co ion is removed, the parameterization for 
the exchange interaction remains the one calculated for 



6 



R\h 


2 


5 


10 


20 


30 


40 


60 


2 


16(25,50) 


18(63,88) 


21(125,150) 


20(251,276) 


20(377,402) 


21(502,527) 


23(754,779) 


5 


26(157,219) 


46(392,314) 


46(785,471) 




44(2355,1099) 






10 


102(628,753) 


122(1570,942) 


135(3140,1256) 


135(6280,1884) 


135(9420,2512) 


138(12560,3140) 


138(18840,4396) 


15 


122(1413,1601) 


152(3532,1884) 


160(7065,2355) 


167(14130,3297) 


170(21195,4239) 


171(28260,5181) 




25 


128(3925,4239) 


158(9812,4710) 


l'^l(19625,5495) 


181(39250,7065) 








30 


130(5652,6029) 


165(14130,6594) 


183(28260,7536) 


191(56520,9420) 









TABLE II: The calculated Tc for a number cylindrical nanoparticles with different radius, R, and length, h. The radius is 
along the vertical axis and the length is along the horizontal one (both in A). In brackets we report respectively the cylinder 
volume (A^) and surface area (A^). 



the perfectly crystalline CoO phas^^^, so that the atom 
removal has only geometrical effects. 

Our results are presented in Fig. |5] where we show the 
specific heat as a function of temperature for different 
Co concentrations, [Co]. In the same figure we also 
report the extracted Tq 's as a function of [Co] . The most 
notable feature is that the peak in C{T) shifts towards 
lower temperatures and gets broader as Co atoms are 
removed. This means that the presence of defects 
reduces the critical temperature, so that lifting locally 
the frustration does not help in improving the magnetic 
properties. We find that Tq decreases linearly with the 
fraction of vacant sites and extrapolates to Tc = for 
[Co]=30%, which is larger than the nearest neighbour 
percolation threshold. This, however, should not be 
surprising since deviation from the linear dependence 
are expected at low [Co]. Interestingly such a linear 
dependance of Tc on the defect concentration has been 
reported previously in both e xperimental and theoretical 
studies for Ni-Cu alloyP^H 

Finally we observe that none of the structures investi- 
gated present magnetizations greater than 0.003 /xb/Co 
above 100 K. We then conclude that the breaking of 
frustration by means of defects appears to have little 
impact on the magnetization and as a consequence 
perforated nanoclusters do not appear to sustain any 
room temperature ferromagnetism. 



III. CONCLUSIONS 

The focus of this paper was to test the feasibility of 
the proposal that uncompensated spins on the surface of 
CoO nanoclusters embedded in ZnO are responsible for 
the measured magnetic properties of ZnO:Co. Previously 
we have demonstrated that bulk wurtzite CoO displays 
a high degree of frustration so that no net magnetization 
is found even at low temperature. Furthermore the 
Tc calculated from the peak in the specific heat is 
substantially lower than room temperature. 

The present study has revealed that finite sized clusters 
are no more promising. Although the frustration is 
lifted at the interface, we found that usually the finite 
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FIG. 5: [Color on line] Specific heat against temperature 
curves for a perforated nanosphere of radius 19 A and various 
Co concentrations, [Co]. We can clearly observe that by 
reducing the Co content the peak in C(T) moves to lower 
temperature and gets broader. In the inset we show Tc as a 
function of [Co] as extracted from the specific heat. 



size reduces the critical temperature for magnetism with 
respect to its bulk value. This happens regardless of the 
shape of the particle and of the presence of Co empty 
sites. For some relatively large spherical and cylindrical 
nanoparticles we have found a marginal enhancement of 
Tc with respect to bulk, for which we speculate on a 
non-monotonic dependence of Tc with particle size. In 
any case the residual magnetizations originating from the 
finite size remain extremely small at any temperature, so 
that finite particles can hardly be considered at the origin 
of the claimed ferromagnetism of ZnO:Co. 

In concluding we wish to remark once again that our 
model is based on a parameterization of the magnetic 
interaction rooted in DFT calculations for bulk CoO. It 
is indeed possible that such an interaction is substantially 
altered at surfaces, so that much larger exchange 
constants may be found either at grain boundaries or 
at free surfaces or at the interfaces with the substrate. 
This can promote residual ferromagnetism. Furthermore 
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our present model does not include intrinsic defects, complexes with Co which may interact over a long rang 
which can both alter the local magnetic coupling or form 
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